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Abstract 


Rapid fluctuations in the emission of solar bursts may have many 
different origins e*g. the acceleration process can have a pulsating 
structure, the propagation of energetic electrons and ions can be 
interrupted from plasma instabilities and finally the electromagnetic 
radiation produced by the interaction of electrostatic and electromagnetic 
waves may have a pulsating behavior in time* In two separate studies 
(Vlahos, Sharma, and Papadopoulos, 1983; Vlahos and Rowland, 1984) we 
analysed the contitions for rapid fluctuations in solar flare driven 
emission* 


1* Time evolution of a beam-plasma system 

A possible interpretation of the solar hard X-ray bursts is that a 
relativly large flux of non-thermal electrons accelerated near the energy 
release site of a solar flare, stream along the magnetic field lines toward 
the chromosphere. The interaction of the non-thermal electrons with the 
ions in the upper chromosphere produces, via collisional Bremsstrahlung, 
photons with energies from 20 Kev to several Mev. The energetic electrons 
are unstable to the growth of plasma waves. In the linear regime , plasma 
waves with phase velocity v^ will grow with a rate 
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where is the bean velocity, « e is the plasna frequency, n^, are the 

beam and ambient densities, is the beam spread around the beam velocity 

and t the angle between the wave vector and the external magnetic field. As 

o 

the beam driven waves reach a threshold value e.g. W../n A T =[v /v.j the 

tn u e e b 

dispersion characteristics of the medium change dramatically and lead to 
the formation of solitons. The long wavelength high frequency waves are 
localised in space, producing soli ton type formations in the wave energy 
and cavities in the ion density. The local ion density is propotional to 
the plasma wave energy (Papadopoulos, 1975) e.g. 
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where is the energy density of the beam exited waves. The localization 

in space of the plasma waves causes a fast transfer of wave energy out of 
resonance with the beam electrons and into resonance with the low energy 
tail of the Maxwellian distribution. The rate with which the waves are 
transfered away from resonance with the beam is 

, nl =(„/M ) 1/2 {w|/n 0 T e ) 1/2 » e (3) 

1 2 

for (W L /Oq^) > ( v e / v jj) • The net result from the soliton formation is: 
(1) the transfer of energy from the beam driven waves to non-resonant 
plasma waves (W^) with low phase velocity and (2) the excitation of 

ion-density fluctuations. As long as (Wg/njjT^ 6n/n Q the ion fluctuations 

are undamped. The non-linear frequency shift, due to W* , for the low 
frequency waves is -kc g , so we obtain an essentially stationary non-linear 

ion pertrubation. When Wg is absorbed by Landau damping on the electrons, 

this non-linear frequency shift disapears, and the ion fluctuations are 
converted to normal linear modes (ion acoustic waves) of the plasma and can 
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then damped. 


The evolution of a beam-plasma system then follows three stages: (1) 
the linear growth of beam driven waves with rate (2) Once the wave 

energy passes a threshold value, discussed above, a fast transfer of plasma 
waves to lower phase velocity waves occurs. (3) Ion density cavities are 

also formed once (wj /n^T^) is above (v^/v^)^* When Wg < W s the 

non-resonant waves cannot support the cavitons and the "cavitons radiate" 
ion density fluctuations. The ion density fluctuations scatter the high 
frequency waves to lower phase velocityies with a rate 

“nl' <V n oV<VvA (4) 

where A =10A_ and A n is the electron Debye length. The time dependent 
s D ue 

evolution of the beam driven waves can also be described by the folowing 
phenomenological system of rate equations 


ar= ’l"! -t nl w 2 > <' , r"2» - “ml"! *<v2> 

(5) 

5T= ’nl* 2 «Cr W th> -V* - a NL W 2* ( V*2 ) 

(6) 

dW , . , 

dt 2 * * S * * 7 NL W 2 '("l-V- -lV<".-»2> 

(7) 


where is the damping of the Langmuir waves in the tail, is the 

damping of the ion density fluctuations and 9 is a step function (see more 
discussion on egs(5)-(7) in Vlahos and Rowland, 1984). We solved eqs(5)-(7) 
numerically and the results are ploted in Fig. (1). 


2. Stochastic three wave interaction 


Electron beams, formed during the impulsive acceleration phase of a 

flare, stream continuously toward the chromosphere. Since the magnetic 

field increases slowly in the direction of propagation a fraction of the 

beam electrons precipitate while the remainder became part of the stably 
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Figure 1: The rate equations are solved numerically using n t/ n o =i0 * 

wj(t=0)=5*10 W2(t=0)=W fl (t=0)=10 v^/v^lO. (a) The temporal 

evolution of the resonant high frequency waves and the acoustic wave 
(b) the temporal evolution of the non-resonant waves. The individual 
pulse have a duration » nsec to msec, dependig on the beam strength. 


trapped component. For a plasma with cyclotron frequency (12 e ) larger than 

the plasma frequency (» e ), the presence of the beam can amplify waves in 

both , lower and upper hubrid brances. The beam exited electrostatic waves 

are convectivly amplified and propagate toward the cutoff region where 

their energy piles up in a narrow wave packet with bandwith 6w trtI . 

UH 

Coherent upconversion of the beam amplified electrostatic waves to 
electromagnetic waves produce an intrinsically stochastic emission 
component. For a coherent interaction we average the basic three wave 
coupling equations over space, retaining however their dynamic character. 
We thus have a zero-dimension but dynamic model. The energy input is steady 
and represents the rate at which electrostatic waves convert into the 
interaction volume. The three wave equations for the real amplitude(a.) and 

phase are 


c - (p 
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dt 1=7 i a r il i a 2 a 3 coat 


4*2 

"dt“ = 1, 2 < *2 _ ' 1 l a 3 a l COst <9) 


da„ 

~dt~ = 7 3 a 3 +/1 2 a l a 2 COSlf> (10) 

dt 

~dt~ = "^UH +( Vl a 2 /a 3" il l a l a 3 /a 2‘ il l a 2 a 3 /a l )sin * (11) 

where are the coupling coefficients and Vj are the linear growth (or 

damping) of the wave. 

The evolution of the electromagnetic wave is shown in Fig. (2). 



t (msec) — " 


Figure 2: The time evolution of the electromagnetic wave resulting from the 
dynamic interaction of two electrostatic waves near the cutoff region. 
The individual pulses have a duration < 1msec. The dotted line 
represents the response of an instrument with resolution >10 msec. 
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